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ABSTRACT 

The Swift-detected CRB 060614 was a unique burst that straddles an imaginary divide between 
long- and short-duration gamma-ray bursts (GRBs), and its physical origin has been heavily debated 
over the years. Recently, a distinct very-soft E814W-band excess at t ^ 13.6 days after the burst was 
identified in a joint-analysis of VLT and HST optical afterglow data of CRB 060614, which has been 
interpreted as evidence for an accompanying Li-Paczynski macronova (also called a kilonova). Under 
the assumption that the afterglow data in the time interval of 1.7 — 3.0 days after the burst are due to 
external forward shock emission, when this assumption is extrapolated to later times it is found that 
there is an excess of flux in several multi-band photometric observations. This component emerges 
at ^4 days after the burst, and it may represent the first time that a multi-epoch/band lightcurve of 
a macronova has been obtained. The macronova associated with CRB 060614 peaked at t < 4 days 
after the burst, which is significantly earlier than that observed for a supernova associated with a long- 
duration CRB. Due to the limited data, no strong evidence for a temperature evolution is found. We 
derive a conservative estimate of the macronova rate of ~ 16.3lg®2^ Gpc“®yr“®, implying a promising 
prospect for detecting the gravitational wave radiation from compact object mergers by upcoming 
Advanced LIGO/VIRGO/KAGRA detectors (i.e., the rate is TZgw ~ Mpc)® yr“®). 

Subject headings: gamma-ray burst: individual (GRB 060614) — radiation mechanisms: thermal — 
binaries: general — stars: neutron 


1. INTRODUCTION 

It is widely accepted that the merger of a binary 
compact object system (either a neutron-star, NS, bi¬ 
nary, or a stellar-mass black hole, BH, and NS binary) 
produces the high-energy y-ray emissio n in a short- 
durat i on gamma-ray burst (SGRB) ev e nt (lEichler et al. I 
119891 : iNaravan et al. I 1199^ iBereer I 1201411 . Indirect 
evidence for SGRBs originating from compact-binaries 
(iGehrels et al. 1 120051 Fong et a l. 1 12010 1 iLeibler et al. I 


120101 : lEong fc Berger I 2013 : IBerger 1 1201^ include the 


location of SGRBs in elliptical galaxies, no associated 
supernova (SN), large galaxy offsets (> 100 kpc) that 
match population synthesis predictions for compact bi¬ 
naries, and weak spatial correlation of SGRBs and re¬ 
gions of star formation within their host galaxies (when 
the hosts can be unambiguously identified). 

A “smoking-gun” signature for the compact-binary 
origin of an SGRB would be the detection of the so- 
called Li-Paczynski macronova (also called a kilonova), 
which is a near-infrared/optical transient powered by 
the radioactive decay of r-process material synthe- 


event (e.g.. Li & Paczvnski 1 

iKulkarni 1 [20f)ft 

Rosswog 1 120051 Metzger et al. 

120101 IKorobkin et al. 

20121 Barnes & Kasen 1 2013 

: iKasen et al. 1 

2013; 

Tanaka & Hotokezaka 

20131: iTanaka et al. 1 

2014; 

Grossman et al. 1 12014 

: iKisaka et al. 1 l20^ 

L5allb: 

LioDuner & Roberts 1 1201,511. 

Macronovae are ex- 


pected to peak in infrared bands, and they display 
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very soft spectra. As such, macronova signals are 
very hard to detect. A breakthrough was made in 
June 2013. In the late afterglow of the canonical 
SGRB 130603B (z = 0.356), an infrared transient was 
interpreted as a macronova produced during a compact- 
binary merger (iTanvir et al. Il201 3t iBerger et al. 1120131 

iHotoicezaka et al. 1120131 iPiran et al. I I20 i4 l Wry re¬ 

cently, a significant F814W-band excess component was 
reported in a re-analysis of the late time optical afterglow 
data (Yang et al. 2015; Y15 hereafter) of the peculiar 
event GRB 060614 that shares some pro perties of both 
long- d uration and sh o rt-dur at ion GRBs (|Gehrels et al. 


20061 iFvnbo et al. I 20061 IPella Valle et al 
Gal- Yam et al. 1 120061 Izhang et al.l 1200'^ The photo- 


;t_aLj 

200ft 


metric spectral energy distribution (SED) of the excess 
component is so soft (the effective temperature is below 
3000 K) that a very-weak SN origin has been strongly 
disfavored. Instead, the excess flux can be interpreted 
as a macronova powered by the merger of a stellar-mass 
BH with an NS (Y15). 

To date, the published literature regarding photo¬ 
metric evidence of macronovae associated with SGRB 
130603B and long-short burst GRB 060614 are based 
on only a single datapoint in each eveniQ. To 
better reveal the physical processes giving rise to 
the macronova emission, multi-band photometric (and 
ideally spectroscopic) observations of the transient 

^ In Fig. 1 of Y15, at f ~ 7.8 days after the burst there was an 
7-band data point that was in excess of the extrapolated power-law 
decline of the hypothesized forward shock emission. However, its 
significance was below 3 (j. 
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are needed in order to com pare observations with 
theoretical predictions (e.g ., [Barnes fc Kasen I 2013 : 
Tanaka fc Hotokezala I [20131: iHotokezaka et al. I 2013; 
Tanaka et al. 112014 iGrossman et al. J12Q141). We there¬ 
fore revisited the extensive data set of GRB 060614 
obtained with the Very Large Telescope (VLT) and 
Hubble Space Telescope (HST) to produce multi-band 
lightcurves and SEDs, and use these information to pro¬ 
vide evidence, and constrain the nature, of the accompa¬ 
nying macronova. 

This work is structured as follows: In Sec. 2 we 
first review the basic assumptions made in Y15. Next, 
we discuss the necessity/feasibility of relaxing these as¬ 
sumptions, and then extract the lightcurve of the associ¬ 
ated macronova. The rate of macronova/compact-object 
mergers is estimated in Sec. 3, and our results and dis¬ 
cussion are presented in Sec. 4. 

2. EXTRACTING THE LIGHTCURVE OF MACRONOVA 
ASSOCIATED WITH GRB 060614 

To robustly establish the existence of a distinct HST 
F814W-band excess in the late afterglow of GRB 060614, 
Y15 assumed that all of the VLT data were due to the 
forward shock (FS) and subsequently fitted the VRI data 
at t > 1.7 day with the same decline rate. In such an ap¬ 
proach, only one F8I4IV-band point at about 13.6 days 
was found to be more than 3a in excess of the fitted 
FS emission. However, the fitted residuals in Fig. 1 
of Y15 display an interesting general trend: the earlier 
data (t < 4 days) were usually negative (with respect to 
the FS afterglow model), while the later data were pos¬ 
itive, indicating that the intrinsic FS emission decline 
was likely steeper than that assumed in their model, and 
there was likely to be an excess of emission at times ear¬ 
lier than 13.6 days. On the other hand, in numerical 
simulations, macronova optical emission usually peaks in 
a few days to a week (rest frame) after the merger event, 
and its subsequent contribution to the afterglow emis¬ 
sion can be non-negligible (e . g.. iBarnes fc Kasen 1 120131 
iTanaka fc Hotokezaka l I2013I iTanaka et al. 11201411 . Af¬ 
ter having solidly established the existence of an excess 
of flux in the analysis performed by Y15, we sought to 
improve the analysis by considering a possible time evo¬ 
lution of the macronova component and modelling the 
entire afterglow dataset accordingly. 

GRB afterglows are expected to be powered by FSs 
that produce synchrotron emission, which have a power- 
law like behavior in both time and frequency, where the 
temporal and energy spectral indices, a and /I, respec¬ 
tively, are defined by fi, oc where t is the time 

since the GRB was first detected by a satellite (e.g., 
IPiran |[200l : iKumar fc Zhang II2015I1 . For SGRBs, the af¬ 
terglow emission emitted after several hours should con¬ 
sist of radiation coming from both the FS and the asso¬ 
ciated macronova. Hence a macronova lightcurve can, in 
principle, be “self-consistently” obtained through a joint 
fit of the observational data. A key outstanding prob¬ 
lem is that current theoretical macronova calculations 
still suffer from significant uncertainties. For example, 
the role of radioactive heating due to the fission of heavy 
r-process nuclei to the energy deposition rate at, for ex¬ 
ample, t ~ 10 days after the merge r , is still poorly under- 
stood fe.g.. iKorobkin et al. Il20ia IWanaio et al. Il20l4 
Hotokezaka 2015 private communication). Moreover, 


the poorly-constrained electron fraction (Ye), the escape 
velocity distribution and the anisotropy of the outflow 
play additional roles in shaping the macronova e mission 
(jTanaka et al. l[2014HLiDDnner fc Roberts l[2015h . are all 
caveats that should be considered when interpreting our 
results. 

In this work we extracted the possible macronova 
emission by decomposing the FS emission from the 
observational data. A reliable estimate of the FS 
emission is very crucial, thus the following facts were 
taken int o account: (i) there was a jet break around 
1.4 days (iDella Valle et al~1l2006l: fMangano et al. Il2007t 


iXu et al.T 2009f) . hence data after this time need only 

be considered; (ii) at t ^ 1.7 — 1.9 days after the burst, 
the optical to X-ray spectrum is well d e scribed by a sin¬ 
gle power-law dPella Valle et al. 1120061 iMangano et al. I 
l2007t iXn et al I l2009fl . suggesting that any macronova 
contribution to the observed flux is negligible; (iii) in 
the interval of 1.7 — 3.0 days after the burst there were 
two measurements in VLT VI bands and three measure¬ 
ments in VLT R band, allowing us to obtain a relatively 
reliable estimate of the FS emission decline. Therefore 
in this work we adopt the VLT and HST observational 
data reduced in Y15, but we assumed that only the VLT 
data in the interval of 1.7 — 3.0 days are due to only FS 
emission, and we used these data to determine the single 
power-law decline of the afterglow. 

The observed magnitudes were first corrected to the 
magnitudes in the R band, assuming an Galactic ex¬ 
tinction of Av=0.07 mag (Schlegel et al. 1998; Schlafly 
& Finkbeiner 2011), the extinction of the host galaxy 
is SMG like Ay=0.05 mag, and the intrinsic afterglow 
spectrum is well described by a single power-law with 
/3 = 0.81 ± 0.0 8, as based on the optic al and UV data at 
150ks fitted bviMangano et al. I (j2007l l and confirmed by 
iXn et al.~l (j2009[ l. The fit to the VLT data collected in 
the time interval of 1.7 — 3.0 days yields a = 2.55 ± 0.09. 
This is steeper than the decay index of a = 2.30 ± 0.03 
obtained in Y15 by assuming all VLT data were FS emis¬ 
sion. Such a difference is reasonable/expected since the 
“underlying” macronova emission contributes to obser¬ 
vations at later epochs, thus causing the LGs to appear 
to decay at a slower rate. In the slow cooling of a jetted 
outflow with significant sideways expansion, when the 
observational frequency is between the so-called typical 
synchrotron radiation frequency t'm and the cooling fre¬ 
quency Uc, the decline and spectr al indexes are exp ected 


to be Of = p ( after the jet break. iSari et al 


(3 = {p — l)/2 (jPiran 112004 iKumar fc Zhang 


1999 1 and 


2015 ). In¬ 


terestingly the observed j5 = 0.81 ± 0.08 and our inferred 
a = 2.55 ±0.09 are in good agreement with the standard 
afterglow model (i.e. they both predict an electron index 
of p ~ 2.6). 

When we subtracted this FS component from the ob¬ 
servational data we found a significant excess in multi¬ 
wavelength bands at t > 3 days, which may constitute 
the first multi-epoch/band lightcurve of a macronova 
ever recorded. The results are shown in Tab. 1 and Fig. 
[U where the errors include the uncertainties of the ob¬ 
served magnitudes and the FS model uncertainties. Al¬ 
though the dataset is still relatively sparse, there is an 
indication that the macronova emission likely peaked at 
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i < 4 days after the merger even10, whic h is consistent 
with current numerical simulations fe.g.. iKasen et al. I 
l2013HTanaka fc Hotokezaka ll201^lTanaka et al. II2014D . 

In comparison, among GRB-associated SNe, SN 2010bh 
had the most rapid rise to maximum bright ness, peaking 
at t = 8.5 ± 1.1 days after GRB 1003160 (ICano et al. I 

l2011l:lBufano et al. il2012l : [Ol ivares et al. Il2012li . which is 

significantly later than that found here. Therefore, in ad¬ 
dition to the remarkably soft spectrum at about 13.6 days 
after the burst noticed by Y15, the rather early peak of 
the excess components found in this work strongly disfa- 
vors an SN interpretation. In the NS-NS merger scenario, 
iMetzger fc Fernandez I ()2014[1 found that some regions of 
the outflow may be Lanthanide-free, and such material 
will become optically thin within a few days after the 
merger and give rise to optical/UV emission lasting for a 
day or so. Such a scenario can, at least in part, explain 
the early peak in the macronova observed here. A suc¬ 
cessful interpretation of the very significant F814W-band 
excess emission at t ^ 13.6 day is, however, a challenge 
for theoretical NS-NS merger models. 

There are five epochs that consist of two or more fil¬ 
ters, which we have combined into five spectral energy 
distributions (SEDs), and are displayed in Fig. [51 An 
excess of flux (i.e., relative to a single power-law spec¬ 
trum) is clearly visible in the three latter epochs. In¬ 
stead, blackbody spectra provide a reasonable fit to the 
observed SEDs (see Fig. [5]). At t ^ 13.6 days after the 
burst, the temperature is estimated to be 2700^7gQ K. 
At other times, the temperatures are poorly constrained 
(i.e., < 4200 K and 3100 —19000 K at 3.86 and 7.83 days 
respectively). Due to their large errors, it is impossible 
to draw any conclusions regarding a possible temperature 
evolution. 

As pointed out in Y15, the progenitor system was 
likely a BH-NS binary, as these types of merger are 
expected to give rise to “bluer”, longer and brighter 
macronova emission than NS-NS mergers due to more 
ejecta mass and a highly-anisotropic distribution of the 
ejecta material fsee lTanaka et al. ll20l4lKvutoku et al. I 
I2015L and the references therein). For GRB 060614, to 
account for the distinct F814W-band excess at t ~ 13.6 
days, a simple estimate based on the generation of the 
macrono va lightcurve in one B H-NS merger model pre¬ 
sented in iTanaka et al. I (I2014D suggests that the ejected 
material from the merger was ~ 0.1 Mq and moved at a 
velocity ~ 0.2c. In this paper, the F814W-band excess is 
just a bit brighter and the parameters of the ejecta are 
likely similar to those in previous estimate. The peak 
emission of VLT//-band (i?-band) excess is as bright as 
^ 24*** mag (^ 25*^ niag), in agreement with the merger 
macronova model (see Fig. |T]). A reliable interpretation 
of the macronova lightcurve, however, requires dedicated 
numerical simulation studies including a proper convo¬ 
lution of the produced complicated macronova spectra 

^ After extracting the possible macronova emission from the 
data, we find that the macronova was always much fainter than 
the afterglow between 1.7 — 3.0 days after the burst and its con¬ 
tribution was smaller than the afterglow uncertainties. Hence our 
assumption that in the time interval of 1.7 — 3.0 days the emission 
is due to just FS is reasonable. However, due to the decline 

of the FS emission and the shallow decay of the macronova, at t>4 
days the contribution of macronova to the total flux can not be 
ignored any longer. 



2 4 6 8 10 20 40 

Time since burst (days) 


Fig. 1.— The observed lightcurves of the macronova associated 
with GRB 060614. Top: The data points are adopted from Y15 
but just the VLT data in the time interval of 1.7 — 3.0 days are 
assumed to arise soley from FS emission, and the solid lines repre¬ 
sent the fit (oc The simultaneous X-ray emission, retr ieved 

from the UK Swift Science Data Centre lEvans et al. Il2009ll . can 
be fitted by plus a constant flux. A const ant X-ray flux of 

(8 ill 4 ) X 10“^® erg cm“^ s“^ was obtained by IMangano et al. I 
11200711 and was interpreted as the emission from a possible Ac¬ 
tive Galactic Nucleus, or it was simply a statistical fluctuation 
because of the low measured flux that was very close to the detec¬ 
tion threshold of Swift XRT. Simultaneous with the very late/weak 
“plateau-like” X-ray emission, the HST F814W-band flux drops as 
^—3.2, Qijt a, possible energy injection mechanism. Bottom: 

Significant excess appears at late times. Note that the data are not 
corrected for any extinction, and only “macronova” emission points 
with a significance above 2cr were kept. The dashed lines, adopted 
from Y15, are macronova model lightcurves generated from nu¬ 
merical simulation for the ejecta from a B H—NS merger, with a 
velocity ~ 0.2c and mass M^j ~ O.IMq, bv lTanaka et al7l II2014I ). 
The green and red lines are in R and I bands, and shadows repre¬ 
sent a possible uncertainties of 0.5 magnitudes (Hotokezaka 2015 
private communication). The macronova model is in agreement 
with the observed data, including those with large uncertainties 
(i.e. significance below 2(t, see Table 1). 

with the response function of the most widely used fa¬ 
cilities in order to aid the comparison with actual data, 
which is beyond the scope of this work. 

3. THE RATE OF THE MAGRONOVA AND 
COMPACT-OBJEGT MERGERS 

So far, two macronovae have been observed at red- 
shifts of 0 = 0.356 a nd z = 0.125 for G RB 130603B 
([Tanvir et al. I 120131 iBerger et al. I I2013D and GRB 
060614 (Y15), respectively. Both events were found to 
be collimated with a half-opening angle 9\ 0.1 (e.g., 

iXu et al.~ll2009t[Fan et al. 1120131 iFong et ai.'ll2014D . For 

the macronovae at z > 0.4, HST observations are cru¬ 
cial to get the signal but HST obser vations of suc h 
“high”-z short GRBs were very rare ([Berger I I2014D . 
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Fig. 2.— The observed SED evolution of GRB 060614. From top 
to bottom are the SEDs at t = (1.73, 2.84, 3.86, 7.83, 13.6) days, 
resp ectively. Solid circl e s are from Y15, red crosses are VLT data 
from l^ella_Valle_et__aL 1 II2006I1 and purple crosses are Swift UVOT 
data from IMangan^^tal. I ||200W) . all have not been corrected for 
extinction. In two early observations the SEDs can be fitted with 
a single power-law spectrum with extinction of the Galaxy and 
the host galaxy, where dash dot line is the intrinsic spectrum and 
dashed lines are extincted. The remaining three observations are 
fitted by a single power-law and a blackbody spectrum (T = 2700 
K, dotted line), where extinction has been taken into account. 


Table 1. The macronova component of GRB 060614 


Time from GRB 
(days) 

Filter 

Magnitude® 

(Vega) 

7.828 

VLT V 

(25.6±0.6) 

3.869 

VLT R 

(25.3±0.6) 

4.844 

VLT R 

24.9±0.3 

6.741 

VLT R 

25.3±0.3 

10.814 

VLT R 

(26.5±0.8) 

14.773 

VLT R 

(27.2±1.0) 

3.858 

VLT I 

23.7±0.4 

7.841 

VLT I 

24.6±0.4 

13.970 

HST F606W 

26.9±0.4 

13.571 

HST F814W 

25.05±0.12 


Note: a. The magnitudes of the extracted macronova component. 
The observations with errors larger than 0.5 mag have been brack¬ 
eted. 

The current sample can be taken as that recordable by 
Swifts an instrument with a field-of-view of 2 steradi- 
ans, in the last ten years for the events with 2 < 0.4. 
With these numbers in mind, we estimated the local 
macronova/compact-object merger rate to be 

7^mac^onova ~ 16.3;l®2" GpC-^ yr-y^j/O.l)-". 

Note that this rate is corrected for beaming, as such 
it is compatible with the un-beamed SGRB rate of 
4 ± 2 Gpc“^ yr“^ (jWanderman fc Piran |[2M^ . For the 
upcoming Advanced LIGO/VIRGO/KAGRA detectors 
that can detect the gravitational wave radiation from 
compact-object rn ergers within a distance D ~ 200 Mpc 
(jAasi et al. Il2013ll . the detection rate is expected to be 

7^GW - Mpc)3 yr-i. 

Bear in mind that such rates are (very) conservative since 
(1) macronova searches usually need HST-like detectors 
and not all SGRBs and long-short GRBs have been fol¬ 


lowed down to deep limits; ( 2 ) it is likely that just a 
fraction of compact-object mergers can produce GRBs. 
Hence the above estimates are better taken as lower lim¬ 
its. We thus conclude that the prospect of detecting 
gravitational wave radiation from merger events in the 
near future is quite promising. 

Interestingly, a realistic e stimate of the BH-NS merger 
rate is ~ 30 Gpc“^yr“^ (lAbadie et al. 1 1201011 . which 
is compatible with T^-macronova estimated here, implying 
that a BH-NS merger origin for GRB 060614 is indeed 
plausible. 

4. DISCUSSION 

Since iLi fc Paczvnski I (|1998ll first proposed that there 
may be a near-infrared/optical transient following the 
merger of a compact binary, significant progress has been 
made in numerical simulations (e.g., [Barnes,^ Case^ 
20131 li^sen et al. 1 120131 iTanaka fc Hotokezaka 1201.31 

Tanaka et al. 112014: iKvutoku et al. II201,3L and the ref¬ 
erences therein). Gonversely, observational macronova 
signatures have only been detected for SGRB 130603B 
(iTanvir et al. ll 2 Ma lBerger et al. II2013I1 and long-short 
burst GRB 060614 (lYang et al. II2015I1 . 

Due to the lack of (detailed) lightcurves and spectra, 
the knowledge we can learn is rather limited and the 
predictions made in the num erical simulations can not 
be fully tested. For example, iHotokezaka et al.~l (|2013[1 
showed that for the single macronova data point of SGRB 
130603B the NS-NS and BH-NS merger scenarios can 
not be distinguished. In the present work, with the as¬ 
sumption that the afterglow data in the time interval of 
1.7—3.0 days after GRB 060614 are generated by external 
forward shock we have shown that at late times there are 
significant excess components in multi-wavelength photo¬ 
metric observations (see Fig. [1]). There is evidence shows 
that the associated macronova likely peaked at t < 4 
days after the 7 — ray transient, which is consistent with 
current numerical simulations of macronova emission but 
much earlier than the peak times of GRB-associated SNe. 

In the approximation of a thermal spectrum, the tem¬ 
perature of the excess component is inferred to be ~ 2700 
K at t ~ 13.6 day. Due to the limited data, no strong 
evidence for evolution of the temperature can be es¬ 
tablished. Such a temperature is significantly lower 
than that of an SN at the same time scale, typic ally 
^ 0.5 — 1 X 10^ K (see e.g., iDella Valle et al. 1120061 and 
iCano et al. 1 1201 ill , but it is similar to that expected at 
the photosphere for the r ecombination of Lanth anides 
(i.e., T ^ 2500 K, see e.g., iBarnes fc Kasen 11201.31 1. This 
lends additional support to the neutron-rich nature of 
macronova ejecta. 

We conservatively estimated the macronova rate 

TT-macronova ~ 16 . 318 ^ 2 ^ Gpc“^yr“\ and implied that 

the detection prospect of the gravitational wave ra¬ 
diation from compact object mergers by upcoming 
Advanced LIGO/VIRGO/KAGRA detectors is quite 
promising, where the expected rate is TT-gw ~ 
0-5lo;L(^/200 Mpc)3 yr-i. 

In the foreseeable future, it is anticipated that in¬ 
creasingly more macronova lightcurves will be recorded. 
There could be two types of macronova lightcurves at 
least. One group is from NS-NS mergers and the other 
is from BH-NS mergers. On one hand, BH-NS merg¬ 
ers are expected to give rise to “bluer”, longer and 
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brigh ter macronova emission than the NS-NS mergers 
ie.g.. iHotokezaka et al. ir2013D . On the other hand, the 
BH-NS merger rate is generally expected t o be at most ^ 
1/10 time that of the NS-NS merger rate (jAbadie et al. I 
1201011 . Hence, if the macronova associated with GRB 
060614 indeed arose from a BH-NS merger, its lightcurve 
will be different from the majority of the sample that is 
expected to be from NS-NS mergers. 
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